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Abstract The present paper reports an inductor-free 
realization of Chua's circuit, which is designed by suit- 
ably cascading a single amplifier biquad based active 
band pass filter with a Chua's diode. The system has 
been mathematically modeled with three-coupled first- 
order autonomous nonlinear differential equations. It 
has been shown through numerical simulations of the 
mathematical model and hardware experiments, that 
the circuit emulates the behaviors of a classical Chua's 
circuit, e.g., fixed point behavior, limit cycle oscilla- 
tion, period doubling cascade, chaotic spiral attractors, 
chaotic double scrolls and boundary crisis. The occur- 
rence of chaotic oscillation has been established through 
experimental power spectrum, and quantified with the 
dynamical measure like Lyapunov exponents. 

Keywords Active Bandpass Filter • Chaotic electronic 
circuits • Chua's circuit • Chaos • Bifurcation 



1 Introduction 

Design of chaotic electronic circuits offers a great chal- 
lenge to the research community for the last three decades 
[T]-[2]. The motivation for designing a chaotic electronic 
circuit comes mainly from two facts: first, one can 'ob- 
serve' chaos, and can also control the dynamics of the 
circuit by simply changing the physically accessible pa- 
rameters of the circuit, e.g., resistor, capacitor, volt- 
age levels, etc.; second, there are multitude of appli- 
cations of chaotic electronic oscillators starting from 
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chaotic electronic communication to cryptography [3]- 
[4 . A Chua's circuit is the first autonomous electronic 
circuit where a chaotic waveform was observed experi- 
mentally, established numerically, and proven theoreti- 
cally [5^-[6^. Moreover, it established that chaos is not 
a mathematical abstraction or numerical artifact, but 
is a very much realizable phenomenon. 

After the advent of chaotic Chua's circuit, a large 
number of works have been reported on different meth- 
ods of realization of this circuit. All of these realizations 
are mainly centered around following goals: inductor- 
free realization of Chua's circuit, and realization of Chua's 
diode. The reason behind the inductor-free realization 
lies in the fact that the presence of an inductor makes 
the circuit bulky, unsuitable for IC design, less robust, 
etc. In the inductor-free realization, the inductor in 
Chua's circuit is replaced by a general impedance con- 
verter (GIC) that requires at least two op- amps. An- 
other approach to this end is Wien-bridge based Chua's 
circuit variant [7] , where a Wien-bridge oscillator is cas- 
caded properly with the Chua's diode. In this context, 
Ref. [9' and fTO^ report a general algorithm for design- 
ing Chua's circuit, in which it has been shown that a 
sinusoidal oscillator can be converted into Chua's cir- 
cuit by incorporating a proper type of nonlinearity. Dif- 
ferent realizations of Chua's diode, which is the only 
locally active nonlinearity in the circuit, have been re- 
ported, e.g., Chua's diode using VGA ^l , [8 , CFOA 
[12 , IC realization of Chua's diode [13], etc. Also, four 
element Chua's circuit has been reported in Ref.[l4], 
which is the minimum component Chua's circuit till 
date. A detailed overview on the Chua's circuit imple- 
mentation can be found in Ref. [15 and [16 . Recently, 
an inductor-free realization of the Chua's circuit based 
on electronic analogy has been reported [8] , [TT] , which 
provided the advantage of exploring the system dynam- 
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ics in the negative parameter space (i.e., negative values 
of Oi and /3 of a classical Chua's circuit [5,). 

Present paper reports an inductor-free realization 
of Chua's circuit, in which we have properly cascaded 
a single amplifier biquad based active band pass fil- 
ter with a Chua's diode. The system has been math- 
ematically modeled with three-coupled first-order au- 
tonomous nonlinear differential equations. It has been 
shown through numerical solutions of the mathemati- 
cal model and real world hardware experiments that the 
circuit shows all the behaviors of a classical Chua's cir- 
cuit, e.g., fixed point, limit cycle formation, period dou- 
bling cascade, chaotic spiral attractors, double scrolls, 
and boundary crisis. Occurrence of chaos has been es- 
tablished through Lyapunov exponents and experimen- 
tal power spectrum. 

The paper is organized in the following manner: 
next section describes the proposed circuit and its math- 
ematical modeling. Numerical simulations, computations 
of nonlinear dynamical measures, e.g., Lyapunov ex- 
ponents are reported in section [3] Section [4] gives an 
account of the experimental results. Finally, section [5] 
concludes the outcome of the whole study. 



2 Proposed circuit and its mathematical 
modeling 

The proposed circuit is shown in Figjlja). This circuit 
has two distinct parts: (i) a second order narrow band 
active bandpass filter (BPF), and (ii) parallel combi- 
nation of a grounded capacitor (C2) and Chua's diode 
{Nn). V\ node of the BPF is connected to the paral- 
lel combination of capacitor C2 and Chua's diode (A^i?) 
through a passive resistor R. Note that, the inductor- 
capacitor parallel combination of a classical Chua's cir- 
cuit has been replaced by a resonator circuit, which is, 
in the present circuit, an active BPF; resonant (or cen- 
ter) frequency of the BPF can be controlled by simply 
varying the resistors R\ or/and (instead of varying 
a capacitor, or an inductor). 

To keep the number of active component lesser we 
have chosen the single amplifier biquad based narrow 
band active BPF proposed by Deliyannis and Friend 
[18] , [19] , which consists of only one single amplifier (in 
the form of an op-amp), two capacitors having same 
values (C), and four resistors (i?i, i?2, Ra^ andi?^). This 
active BPF is a second order system with the following 
transfer function: 



H{s) 



-{k + l)s/RiC 



s2 + {2/R2C - k/RiC)s + (l/i?ii?2C2) ■ 



(1) 



Here, k = R^/Ra. Proper choice of Ra and R^ (and 
hence k) makes the coefficient of s in the denominator 



negative, which in turn brings the pole of the circuit 
to the right half of the s-plane, resulting a sinusoidal 
oscillation. Real time dynamics of the system can be 
expressed in terms of the following two coupled first- 
order autonomous differential equations ^20] : 



C 



dVi 



Ri 



(2fe + l) 
" {k + l)i?2 

at Ri R2 

Equation ^ can be written as: 
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It can be seen from (|3| that the BPF can be made 
to oscillate sinusoidally if one can ensure the follow- 
ing condition: an + 0^22 = 0. This in turn gives the 
condition of oscillation of the sinusoidal oscillator as: 
k = 2R1/R2; subsequently, the frequency of oscillation 
of the sinusoidal oscillator can be obtained from the 
condition:y'aiia22 — ci^i2<^2i = 0. Thus the frequency 
of oscillation is obtained as: uq = l/C^/^i^- 

Chua's diode is characterized by the function /(V2), 
where V2 is the voltage drop across C2. /(^2) is a three 
segment piece wise linear function with Ga and G5 as 
the slopes, and Bp is the break point voltage for those 
segments. Variation of f{V2) with V2 has been shown in 
Figjljb) (showing only the three segment region). /(V2) 
is defined by: 



f{V2) = GbV2 + 2 (Ga - Gb)i\V2 + Bp\ - \V2 - Bp\ 



(4) 



In the present work, we have chosen the VOA imple- 
mentation of Chua's diode [6 as shown in Figjljc). 

The dynamics of the proposed circuit can be de- 
scribed by three-coupled first-order autonomous non- 
linear differential equations in terms of V2, Vi and Vb 
(Fig.[l]): 



C: 
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dt 
dVi 
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dt 
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R 

(fc + 1) 
R 



- f(V2), 
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(5a) 
(5b) 
(5c) 



Where, i?' = RiR/{Ri + R). 

Eq.([5]) has been written in the following dimensionless 
form using dimensionless quantities: r = t/RC, u = 
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Active Band pass Filt&r 




Fig. 1 (a) Active bandpass filter based Chua's circuit, (b) 
Characteristics of the Chua's diode, (c) VOA implementation 
of Chua's diode iS], [H]. 



{u x,y,z), X V2/Bp, y Vi/Bp, z = Vo/Bp, n 
R/Ru r2 = R/R2.a = C/C2: 

X = a[y — h{x)]^ (6a) 

y = -kx + k{ri + l)y - -^^^^t^y^^, (6b) 

(A: + 1) 

z = -{k^ l)x + (A: + l)(ri + 1)^ - 2r2Z. (6c) 
Here, h{x) is defined as 

h{x) =x + f{x) 

=mix + 0.5(mo — mi){\x -\- 1\ — \x — 1|), (7) 

where, we have defined the following parameters: mo = 
{RGa + 1), and mi = {RGb + 1). 



3 Numerical studies: Phase plane plots, 
Bifurcation diagrams and Lyapunov exponents 

Numerical integration has been carried out on (|6| using 
fourth-order Runge-Kutta algorithm with step size h = 
0.001. We have chosen the resistor Ri as the control 
parameter (remembering that Ri controls the gain and 
center frequency of the BPF); thus in the numerical 
simulations, ri acts as the control parameter keeping 
other parameters fixed at: r2 = 0.2, k = 0.04, a = 20, 
mo = —1/7, and mi = 2/7. Note that, the values of mo 
and mi are generally used in a classical Chua's circuit. 

It has been observed that, with increasing ri, for 
ri < 15.08, the circuit shows a fixed point. For ri > 
15.08, the fixed point loses its stability through Hopf 
bifurcation, and a stable limit cycle emerges. At ri = 
15.54, limit cycle of period- 1 becomes unstable and a 
period-2 (P2) cycle appears. Further period doubling 
occurs at n = 15.72 (P2 to P4), and n = 15.76 (P4 to 
P8). Through period doubling bifurcation, the circuit 
enters into the regime of spiral chaos at ri = 15.78. 
With further increase of ri, at ri = 17.02, the cir- 
cuit shows the emergence of a double scroll attr actor. 
Finally, the system-equations show diverging behavior 
beyond ri = 24.81, indicating boundary crises. Phase 
plane representation {in y — x plane) for different ri is 
shown in Fig|2j which shows the following characteris- 
tics: period-1 (ri = 15.5), period-2 (ri = 15.6), period-4 
(ri = 15.75), period-8 (ri = 15.77), spiral chaos(ri = 
16.35), double scroll (ri = 17.59). The phase plane 
plots, for the same parameters, in z — y plane is shown 
in Figj3]^a)-(e) also, the same double scroll chaotic at- 
tractor in z — x plane is shown in Figjsjf). 

These observations can be summarized through bi- 
furcation diagram with ri as the control parameter. 
Bifurcation diagrams are obtained by using Poincare 
section ^ at 7/ = 0.1 with ^ < 0, excluding the 
transients. FigjI] (upper and middle trace) shows bifur- 
cation diagrams in x and respectively, for different 
ri. Clearly, it shows a period doubling route to chaos. 
Further, it shows the presence of a period-3 window at 
ri = 15.93. Among the other periodic windows, more 
prominent is a period-5 window near ri = 18 inter- 
spersed in a double scroll chaos. For quantitative mea- 
sure of the chaos generated by the circuit model, Lya- 
punov exponents have been computed using the algo- 
rithm proposed in |22j. Fig. |4] (lower trace) shows the 
spectrum of largest Lyapunov exponent in ri parameter 
space. It agrees with the bifurcation diagram. Presence 
of a positive LE and strange attractors (as the dou- 
ble scroll attractor) ensures the occurrence of chaotic 
behavior in the circuit. 
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Fig. 2 Phase plane representation (in y — x plane) for dif- 
ferent ri: (a)ri = 15.5 (period-1), (b)ri = 15.6 (period- 
2), (c)ri = 15.75 (period-4), (d) n = 15.77 (period-8), 
(e)ri = 16.35(spiral chaos), (f)ri = 17.59 (double scroll). 
(r2 0.2, k 0.04, a 20, mo -1/7, and mi 2/7) 
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Fig. 3 Phase plane plots in z — y plane (a)-(e), and z — x 
plane (f) for different values of ri. Parameter values are same 
as Fig[2l 

4 Experimental results 

The proposed circuit has been designed in hardware 
level on a breadboard. Since the circuit needs three 
VOAs we have chosen IC TL074 (quad JFET op-amp) 
for the present purpose with ±12 volt power supply. 
BPF has been constructed with C = 100 nF, Ra = I 
kQ^ Ri) = 56^1, and R2 = lOkQ. Chua's diode is con- 
structed with following parameters JT]: Rs = 2.2 

kn, i?4 = 220 = 220 Rq = 3.3 kr^, R7 = 22 kn, 

and Rg = 22k^}. The grounded capacitor value is taken 
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Fig. 4 Bifurcation diagram of x (upper trace) and z (middle 
trace) with ri as a control parameter. Lower trace shows the 
Largest Lyapunov exponent (Xmax) with ri. (r2 = 0.2, k = 
0.04, a 20, mo -1/7, and mi 2/7). 

as C2 = 5 nF. Coupling resistor R is fixed at l.Skfl (ap- 
prox.) using a 2 kQ POT. To explore the dynamics of 
the circuit we have varied the resistor Ri through a 
1 kQ POT. All the potentiometers are precision POT 
having thousand turns. Capacitors and resistors have 
5% tolerances. 

For Ri > 200Q the circuit shows a fixed dc value 
(equilibrium point of the circuit). For Ri < 200Q a 
stable limit cycle has been observed with frequency 
1955Hz. At Ri = 185^1 (approx.) the limit cycle of 
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period- 1 loses its stability and a period-2 oscillation 
emerges. A period-4 behavior has been observed at Ri = 
177Q (approx.). period-8 is found for Ri = 1751^ (ap- 
prox.). Further decrease in Ri results in spiral chaos in 
the circuit (at Ri = 170r^(approx.)). The double scroll 
attractor is observed at Ri = 122Q (approx.). The cir- 
cuit shows a large limit cycle for Ri < 6SQ that indi- 
cates the occurrence of boundary crises, where the ac- 
tive resistor becomes eventually passive , i.e., for a large 
voltage across its terminals, the i — v characteristic of 
Chua's diode is no longer a three-segment curve but, it 
becomes a five-segment curve. Two additional segments 
(situated in two outer regions) give positive resistance 
behavior, i.e., now the instantaneous power consump- 
tion becomes positive All the above mentioned be- 
haviors (except the large limit cycle) have been shown 
in Figjs] (in V1-V2 space) and Figj6] (in Vq-Vi and V0-V2 
space), which depict the experimental phase plane plots 
recorded in a real time oscilloscope (Aplab make, two 
channel, 60MHz). 

Apart from the experimental phase plane plots, an- 
other very useful tool for exploring the real circuit be- 
havior for a range of circuit parameters is experimental 
bifurcation diagram [23], [24]. To grab the behavior of 
the proposed circuit for a range of values of ex- 
perimental bifurcation diagram of Vq has been plotted 
taking Ri as the control parameter (Figj?]). For seventy 
different values of Ri in the range of 195 Q to 120 
we have acquired the experimental time series data of 
Vo using Agilent make Infinium digital storage oscil- 
loscope (5000 data per set) and find out the poincare 
section using the local minima of the time series data. 
Plotting these for all the Ris gives the experimental 
bifurcation diagram. Schematic difference between the 
numerical bifurcation diagram (Figj4] (middle trace)) 
and experimental bifurcation diagram occurs due to the 
difference in the sampling schemes employed in finding 
out the Poincare section. However, qualitatively they 
agree with each other, e.g., both the diagrams show pe- 
riod doubling route to chaos. Blurring of points in the 
experimental bifurcation diagram occurs due to the in- 
herent circuit noise that is reflected in the experimental 
time series data. Also, since we have to acquire the data 
manually, resolution of the plot is not high enough in 
comparison with the numerical bifurcation plot. 

Time series data and the corresponding FFTs of V2 
are measured using Agilent make Infinium digital stor- 
age oscilloscope (500MHz, sampling rate IGs/s) and are 
shown in FigJsJ Figjsj^a) (upper trace is for time series 
and the lower trace is for FFT) shows the period- 1 os- 
cillation with frequency 1955Hz; Figjsj^b) and (c) shows 
the results for spiral chaos and double scroll chaos, re- 
spectively. Power spectrum of Figjsj^b) and (c) (lower 
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Fig. 5 The oscilloscope trace of experimentally obtained 
phase-plane plots in the V1-V2 space (a) Period- 1 at Ri = 
195n (b) Period-2 at Ri = ISOQ (c) Period-4 at Ri = 176^ 
(d) Period-8 at Ri = 175f2 (e) Spiral chaos at Ri = 150n. 
(f) Double scroll attractor at Ri = 120Q. (Other parameters 
are: BPF:C = 100 nF, Ra = 1 kQ, R^ = 56^, and R2 = 10k 
n. Chua's diode: R3 = 2.2 kQ, R4 = 220 R5 = 220 
Re = 3.3 k^^, R7 = 22 k^^, and Rg = 22 k^^. The grounded 
capacitor C2 = ^ nF. R = l.SkQ). (a)-(e)Vi {x -axis): 0.1 
v/div, V2 (2/-axis): 1 v/div. (f)Vi (x -axis): 0.3 v/div, V2 (y- 
axis): 2 v/div. 

traces) are continuous and broad in nature, indicating 
the occurrence of chaotic oscillations. It can be seen 
that the experimental results agree with numerical sim- 
ulations of the mathematical model of the circuit. 



5 Conclusion 

In this paper we have reported a single amplifier biquad 
based inductor- free implementation of Chua's circuit, 
where we have suitably cascaded an active band pass 
filter with a Chua's diode. The proposed circuit has 
been modeled mathematically by three-coupled first- 
order autonomous nonlinear differential equations. Nu- 
merical simulations of the mathematical model of the 
circuit and the real world hardware experiments con- 
firm that the proposed circuit shows the same behavior 
that is shown by a classical Chua's circuit (e.g., fixed 
point behavior, limit cycle formation through Hopf bi- 
furcation, period doubling, spiral chaotic attractor, dou- 
ble scroll attractor, and boundary crisis). The presence 
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Fig. 6 The oscilloscope trace of experimentally obtained 
phase-plane plots in the (a)-(e)Vo-^i space ((a)-(d)Vb (x - 
axis): 1 v/div, Vi (^/-axis): 0.1 v/div. (e)Vb {x -axis): 2 v/div, 
Vi (y-axis): 0.3 v/div.), and (f)Vb-V2 space (Vb (a^ -axis): 2 
v/div, ^2 (i/-axis): 2 v/div.); Parameter values are same as 
used in Figjs] 
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Fig. 7 Experimental bifurcation diagram of Vq (v) with Ri 
(Q) as the control parameter. Other parameters are same as 
Fig|5] 



of chaos has been established through the Lyapunov 
exponents and experimental power spectrum. 

As the circuit is inductor-free it has got all the ad- 
vantages of an inductor-free circuit (e.g. suitable for IC 
design, robustness, etc.). Further, instead of varying an 
inductor or a capacitor, we have used a single resis- 
tor (e.g., Ri) as a control parameter to observe all the 
complex behaviors of the circuit. Since the circuit has 
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Fig. 8 Experimental time series and FFT of V2 for (a) 
Period- 1 oscillation (frequency 1955Hz) (b) spiral chaotic os- 
cillation and (c) double scroll chaos (circuit parameters are 
same as described in Fig[5|. (Frequency span of the FFT: 25 
kHz). 



a large number of choices of control parameters (viz. 
i?2, R and a), one can use these parameters to observe 
different behaviors of the circuit. But, the basic bifur- 
cation routes to chaos remains the same for each of the 
parameters. The BPF circuit we have used is narrow 
band (unlike Wien-bridge based circuit [7], where the 
frequency selective network is not narrow band), thus 
one has more control over the center frequency of the 
circuit. The present circuit can be suitably designed 
(with high frequency op-amps) to generate chaotic os- 
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dilations in a high frequency region, and also may be 
useful for chaos based communication systems. 
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